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Abstract 

The calculation of the spreading width of a compound nuclear reaction caused 

by a symmetry breaking i^'-body force acting in an ^-body system (K <^ A, 

usually K = 2) involves the determination of the local average square matrix 

element in A-hody space. This problem is reduced to finding the global mean 

square matrix element in ET-body space. The result is a compact formula 

for the spreading width which contains as an input. Our method is based 

on the dilute gas approximation for excitons close to the Fermi edge. The 

relative strength of the contributions of operators with different exciton rank 

as well as the connection between the energy dependence of the spreading 

width and the body rank of the underlying interaction are established. 
24.60.Dr, 24.80.Dc, 25.70.Gh 
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I. INTRODUCTION 



The compound nucleus (CN) has received considerable attention — both theoretical and 
experimental — in the recent years because of its seeming ability to enhance the effects of the 
weak interaction to a few percent (see for example [|I| and references therein). In the present 
study we consider the breaking of symmetries caused by i^'-body forces in the compound 
process in general. The adequate quantity to characterize symmetry nonconservation in 
many-body systems is the spreading width due to the underlying interaction [^. Although 
in most applications a two-body force is used, we treat the case of arbitrary K since it will 
allow for some insight into the energy dependence of the spreading width. This provides a 
generalization of the results presented in and the communication of some details omitted 
therein. 

In the statistical theory of CN reactions the matrix elements of the interaction are as- 
sumed to show the characteristics of the Gaussian orthogonal ensemble (GOE). The crucial 
point is that statistical assumptions of this kind can only be made about the defining matrix 
elements, i.e. the matrix elements in i^-body space. However, the quantity of physical in- 
terest is the matrix element in A-body space. The connection between the properties of the 
i^-body matrix elements and those in A-body space is called the propagation of the defin- 
ing matrix elements In the present article a solution for this problem is offered which 
consists of the transition to the exciton picture (with the accompanying simplification of 
the basis states and complication in the description of the interaction), the propagation into 
the subspace of fixed exciton number, and the averaging over subspaces, which implies the 
return to the body picture. This procedure is necessary because our formalism makes use 
of the dilute gas approximation (DGA) which — in contrast to the body picture — is very 
good in the exciton representation of the system. In the final sections then, the general 
expression, which involves a convolution of partial and total densities of states, is evaluated 
by inserting well known analytical formulae for the level densities, the results are discussed, 
and limitations of our approach are indicated. 
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II. CONCEPTS 



A. Spaces and bases 

The physical quantity we are interested in is the spreading width 

r^(E) = 27r<V2>p(E) . (2.1) 

Here, ^ ^ is the mean square matrix element of the -f^-body interaction V, and p is 
the level density of the system. The matrix elements are calculated in a basis of eigenstates 
to those parts of the Hamiltonian that dominate the behaviour of the system. The sprea- 
ding width due to additional, symmetry breaking interactions then measures the extent of 
symmetry breaking in the system. This may include the breaking of the independent particle 
structure, isospin symmetry or parity by the residual strong, the electromagnetic or weak 
interaction, respectively. In order to properly treat the variation of with energy, the 
average <^ ^ is limited to states in the neighbourhood of some given excitation energy 
E. In fact, part of the present work will consist of the calculation of the strength function 

S{E\ E) = < > p{E')p{E) (2.2) 

which implies the average over squared matrix elements between configurations close to E 
and configurations close to E' . Once the basis has been specified in detail, this local average 
will be defined in section |ITT| . The basis we work with is built up from the single particle 
states I A) that satisfy the canonical Hartree-Fock equations 

(u + uhf)|A) =: ho|A) = £a|A) , a = 1 . . . D . (2.3) 

Here, D is the dimension of the one-body space spanned by the discrete set of bound 
states, u is the one-body kinetic energy operator, and uhf is the Hartree-Fock mean field 
operator constructed from the strong nucleon-nucleon interaction, which does not include 
the symmetry breaking interaction V for which the spreading width shall be calculated. Let 
us introduce the creation operators and the physical vacuum | ) so that 
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|A) = b[\). (2.4) 

These operators fulfill the fundamental anticommutation relation 

{bx,hi,} = 6xx' . (2.5) 

A basis of A-hodj states is then given by all possible A-fold applications of creation opera- 
tors: 

|A)^|Ai...AA) = bl^...bij ) , (2.6) 

with the condition 

Ai < A2 < . . . < Aa. (2.7) 
The energy of this many-body state is 

^A = EeA, . (2.8) 

i=l 

The energy of the ground state \1 . . . A) is J2t=i and will be denoted £a- The states ( [2.6[ ) 
span the whole A-body space of dimension but for low energy considerations it is 

convenient to split the A-body space into subspaces of increasing complexity and energy 
and to restrict the description of the system to the subspaces in the energy range of interest. 
Technically this idea is realized by the transition to the exciton picture. For moderate 
excitation energies the state of the system does not strongly differ from the ground state 
and there will be only a few states occupied above and unoccupied below the Fermi energy 
Ea- This will make it possible — in section |III A| — to introduce the so-called dilute gas 



approximation and treat the present problem in closed form. Let us therefore introduce the 
exciton creation operators [^-|^ according to ^ 



= bl, for A > C 



a{ = bx for A < C (2.9) 

acting on the exciton vacuum 



|0) = |1...C). (2.10) 

One then has aA|0) = for all A and the anticommutation relation (|2.5|) holds for the 
operators a, too. Briefly: the excitons are fermions. The single body space is split into two 
subspaces, the single particle space with dimension dp = D — C and the single hole space 
with dimension = C. The exciton vacuum |0) — the core — need not be identical to the 
ground state of the A-body system under consideration, i.e. at present we do not fix 

A = A-C. (2.11) 

The energy of the vacuum is called 

c 

Sc = J2e^. (2.12) 

A=l 

Every configuration |A) can then be characterized by a vector 

P={n,...,rp) , r, >C (2.13) 

of its particle configurations and by a vector 

H = {ai,...,ah) , ai<C (2.14) 

of its hole configurations. Every independent particle configuration |A) can be expressed in 
the form 

|^)^at^...at^at^...atjO). (2.15) 
The combination of operators appearing in this equation will also be written as 

at^...at^at^...at,^ = (A^)t. (2.16) 
As usual we introduce the energies 

er,=er,-€c (2.17) 

of the particle configurations relative to the core level as well as the energies of the hole 
configurations 



ea, =ec- Sa, . (2.18) 
The energy of the state ( |2.15| ) then is 

h p 

£P = Y.^a,+T. +^c + ec^ = + £c + ^cA . (2.19) 

i=l i=l 

The symbol has been introduced to simphfy the notation. In the present paper the letters 
r, s, t,u or V will be used for particle states and a, (3, 7, 5 or e for hole states; configurations 
of the physical constituents of the system — called bodies — are labelled A, /x, u, p or a. If 
one chooses C < A, the ground state of the system has A particles and no holes, otherwise 
it has no particles and —A holes. Generally one has 

p-h = A (2.20) 

and 

P > Pmin = max(A, 0) , h> max(-A, 0) . (2.21) 

The maximum number of particles is A and that of holes C = A — A. Hence, the exciton 
picture provides a decomposition of the A-hodj space into mutually orthogonal subspaces Up, 
P = Pmin, ■ ■ ■ ,A. The dimension of the subspace Up is d{Up) = (^'^1^^ (^'^^ . We use the following 
notation for many-exciton states: A state | contains the same numbers of particles and 
holes as | ^) on possibly different single exciton states. A ket | ^) on the other hand differs 
from I ^) with respect to exciton number as well as single exciton states. 



B. Conventions of summation 



The spectral average over the squared matrix elements of V implies sums over the basis 
states ( |2.15|) . To avoid double counting of the basis states a definite order of the indices r 
and a in eq. ( p.l5|) must be observed. We therefore introduce the notation 



E= E • (2-22) 

{H} ai<...<«h 



Here, the a, with i = 1, . . . ,h run over the configurations 1 < < C. The sum is called 
restricted because of the restriction ai < . . . < ah imposed on the permitted terms. Hence, 
the sum has (^"^^^ terms. In contrast the sum 

E= E (2-23) 
is called unrestricted and has (dh)^ terms. Corresponding conventions are used for the set 



P of particle configurations. For later purposes we compress the notation of (|2.22|) further: 



E -E---E' (2-24) 

{HP;...;LT} {HP} {LT} 

which means that the ordering has to be observed only within each group of indices separated 
by semicolons. In the following use will be made of the identity 

E(I|V|^)^ = ^E(I|V|S)% (2.25) 

{H} H 

which holds because the squared matrix element is completely symmetric with respect to 
the h indices in H, but vanishes if any two of them coincide. 



C. Interactions for bodies and excitons 

Now that the basis we are going to work with is specified, we turn to those parts of 
the interaction that have not been considered in its determination and consequently cause 
transitions between basis states. Apart from the residual strong interaction this includes 
the electromagnetic, weak, and other possible forces. In occupation number formalism, a 
i^-body operator 

V= E v(x.,...x.,) (2.26) 
n<...<ix=l 

has the form 

V = E {f^i...fiK\yWi... i^kH^ . . . bt^b,, . . . , (2.27) 
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with the totally antisymmetric matrix element (/ii . . . /ii^|v|z/i . . .vk)- How is this represen- 
tation affected by the transition to the exciton picture? The range of summation of every 
index is split into two parts, body operators are replaced by exciton operators according to 
eq. (|2.9| ), and the resulting terms are brought into normal order and grouped according to 
particle-hole structure. For K = 2, the result is given in refs. and is reproduced in 

Tab. I. The interaction is the sum of the fourteen terms k^{a,q) listed there together with 
their Feynman diagrams The diagrams facilitate the visualization of the systematics in 
and the generalization of the contents of Tab. I, see Fig. |l| . The exciton operators can be 
classified by the three numbers /c, a, and q. The rank k of k^{a,q) is half the number of 
external lines in the corresponding diagram. One finds Q < k < K . The parameter a is the 
number of particle- hole pairs created by k^{a,q) '■ the interaction does not conserve the num- 
ber of excitons; the conservation of the number of physical bodies, however, requires that 
the exciton number changes by particle-hole pairs. The range of a obviously is —k < a < k. 
Finally, q is the number of upward arrows in the diagram. One recognizes that q changes 
in steps of two, since for fixed a every additional particle before the interaction leads to an 
additional particle after the interaction. The range of q is found to be \a\ < q < 2k — \a\. 
Hence, 

K k 2k-\a\ 

'^ = 11Y. (2.28) 

k=0 a=—k <3=kl 

On the table, the operators with k < 2 have contracted hole lines that represent the inter- 
action of the excitons with the nuclear core. It is easily seen (and holds for arbitrary K as 
well) that the number of particle and hole lines before and after the interaction is given by 
N_ = k + a and = k — a, respectively. Out of these there are p± = particles and 
h± = N± — p± = k — holes. We give a few identities that clarify the significance of these 
quantities: 

+ = 2k 
P++P- = (l 
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+ = 2k — q 

p^-p_ = h+-h^ = i(iV+ - iV_) = a (2.29) 
— /i+ = p_ — /i_ = g — A; 
+ = + h_ = k . 

Furthermore, we introduce he = K — k, which is the number of contracted hole indices of 
k^(a,q)- The contracted indices never appear with exciton operators. Since as a consequence 
of ( p.9|) hole indices associated with creators (annihilators) appear in the bra (ket) of the 
matrix element, the general structure of k^{a,q) is ||T0|| : 

fcV(a,g) = J2 J2 J2 J2 J2 

{ai...ahp} {/3i.../3h_} {si...Sp_} {7i---7h+} {ti-tp^} 

(ai . . . ahcf3i . . . (3h_ti ■ ■ ■ tp+ |v|ai . . . ahcli ■ ■ ■'Ih+si . ■ . Sp_) 
X al^^ . . . ala\^ . . . ai,^^a^,_ . . . a^^a,, . . . a,^_ 

Y: {HcH^P+\v\HcH+P.) {A^^X)\a^h-_) . (2.30) 

{Hc;H_P^;H+P+} 



D. Partial and total level densities 

In the statistical model of CN processes, reaction rates are dominated by the available 
phase space. This calls for a detailed knowledge of the densities of states, which have been 
the subject of intensive studies, see for example |TT|-p^. In the present section, we quote 



the results relevant for the sequel. The density of states for the A-body system, 

p{S) = J26{S-S^) , (2.31) 

is approximated by continuous expressions derived with methods of statistical mechanics. 
The famous Bethe formula is 

, exp[2(7r2oE/6)i/2i 
p(B)(^) = yi K y / J — i_ (^2.32) 
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Here, E is the excitation energy of the system, E = S — Sa (with the ground state energy 
Sa) and g is the single body density at the Fermi edge. Blatt and Weisskopf |T5[ give the 
expression 

p(BW)^^) = Cexp[2{aEy/^] , (2.33) 

which may be understood as an approximation to eq. ( |2.32| ) in so far as it contains only the 
term varying most rapidly with energy. Finally the approximation of constant temperature 
(CTA) yields a purely exponential increase for the nuclear level density, i.e. 

P^^^^Hy) = P(yo)e(^-^«)/^ . (2.34) 

Here we have called the excitation energy y, because we will need the total level density in 
this form in section |1V lj|. The nuclear temperature is defined as 



T=J^. (2.35) 
V ^9 



Gilbert and Cameron |]T2[ and v. Egidy and collaborators have compared these 
expressions with experimental data on nuclear level densities and found that the CTA gives 
a good fit up to excitation energies of approximately 10 MeV, whereas above this value the 
Bethe or Blatt and Weisskopf expressions must be used. 

Since we are going to work with the exciton picture, we need expressions for the densities 
of states with fixed exciton number. They will be characterized by the number of holes and 
particles that occur: 

Y: 6iE - E^) ^ p,^^iE) = p,_^^,{E) . (2.36) 

{HP} 

Here, Ef^ = £^ — Sa = ~^ + — £a = ^ ^ ■ The shift S is nonzero only iiC ^ A. 
Summation over the subspaces yields the total density of states: 

E Pp-A,p(^) = P{E) . (2.37) 

p=max(A,0) 

The identity ( |2.37] ) is of course independent of A. A different A merely implies a different 
description, but does not alter the physics of the level density of the A-body system: 
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X: E S{E-E^^^)= E E S{E - E^_^,) = p{E) . (2.38) 

p=max(A,0) {P-A,P} p=max(A',0) {P-A',P} 



Since we are going to make use of this invariance property of p in section |III B| , it is necessary 



to explain in detail how it is to be understood. The sum over p on the r.h.s. of eq. (|2.38|) 



may contain subspaces (e.g. Ip 3h configurations) that do not appear at all on the left 
hand side (which could start with the subspace of 2p Oh configurations). And also the 
single exciton spaces are different: an increase of A decreases the dimension of the single 
hole space and enlarges that of the single particle space. If we assume a constant spacing 
of the single body levels, however, the single exciton energies and in eqs. ( |2.17| ) 



and ( [^.18D that occur in the summations over the subspace configurations are not affected 



by the transformation (|2.38D — except for the highest excited states. We can safely ignore 



this difference in the energy range of interest. Explicitly, the identity (|2.38|) then reads 



E E S{E-e^p_^-S)= E E S{E - e^p^^, - S') , (2.39) 

p=max(A,0) {P-A,P} p=max(A',0) {P-A',P} 



where S' = £c' +^c'^' — ^a-, see below eq. ( |2.36| ). Ericson gave the well-known approx- 
imate expression 

P^^^ - h\p\{h + p-l)\ ' ^^-^^^ 

for the partial density of states ( p.36| ). It is valid if the Pauli principle for excitons is ignored. 
This corresponds to the "dilute gas approximation" discussed below. Eq. (|2.4CI| ) was lateron 
improved in many respects (e.g. [0,|l3]), but is widely used because of its simplicity. 



III. PROCEDURE 



A. Propagation into the Subspace 



The strength function for transitions between states of energy E and exciton number 
N = p + h and states of energy E' and exciton number A^' = t + / is defined as 
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S^t,,){E\E)= Y: CMmE'-El)S{E-E^). (3.1) 

{LT;HP} 

The local average < > over the squared matrix elements between states with E, N and 
E', N' then is 

< > p, ,iE')p,^^iE) = S^t,p)iE', E) , (3.2) 

where the partial densities of states are defined in eq. ( |2.36| ). Hence, the average is defined 
with the delta functions of the energies as weighting factors. However, in the present context 
5{E — E^) is to be understood not as the Dirac distribution but rather as a peaked function 
of suitable width. The "suitable width" is large compared to the average level spacing and 
small compared to intervals over which secular variations of the level densities occur. This 
guarantees that densities of states and strength functions are smooth and can be reasonably 
approximated by the expressions of subsection |IID| . The strength produced by the interac- 
tion operator k^{a,q) is by specialization of eq. (|3.1|) and use of eq. ( |2.30|) : 



^:^{E',E)^ J: {l\ ,yM\H?5{E'-El)5{E-E^^) (3.3) 



{LT;HP} 



X! \L 

{LT; HP} 



Y: {HcH^p^\v\HcH^p^) (a^;)^ (a^:) 

{Hc;H+P+;H_P_} 



^^mE'-El)6{E-E^). 



Note that of the seven quantities that specify the strength only six are independent since 
t = p+a. Furthermore we observe that a nonzero contribution comes only from the subspaces 
with p > p- and h > h^, because otherwise the matrix element vanishes. We proceed to 
estimate this expression in three steps. 

(i) In the first step the i^-body matrix elements of V are considered to be entries of a 
random matrix. Invoking time reversal invariance we postulate that V (in i^'-body space) 
belongs to the Gaussian Orthogonal Ensemble (GOE). This means that the second moments 
of the matrix elements in eq. ( p.27| ) can be expressed by a single parameter v"^, namely (a 
bar over a symbol denotes the ensemble average) 



{KMK2){K[m',) = v\bff^l + (3.4) 

12 



for i^-body configurations i^i, K2. Here, S§, is a generalized Kronecker symbol |T9 

with the properties 

1 : if i^' is an even permutation of K 
—1 : if i^' is an odd permutation of K 



(>K' 



(3.5) 

: if two indices out of K or out of K' coincide 



: if i^'' is not a permutation of K . 
This allows us to estimate the strength function ( |3.3D by its ensemble average, or — by the 
same token — the spectral average < > is identified with the ensemble average. 

(ii) How does the correlation rule (|3.4|) translate into the exciton picture? We find the 
following relation: 



(3.6) 



In principle, eq. ( |3.6| ) allows the simplification of the various terms in the strength func- 
tion ( |3.1| ). The evaluation of the sum over the hole indices {He', H^} in eq. ( |3.3| ), how- 
ever, is complicated by the fact that it is not completely restricted in the sense of eq. ( p. 221 ). 
We therefore use an approximation to eq. (|3.6|): 



{HcH.P+\y\HcH+P^){LcL^T^\y\LcL^T^) = v'S^^Jd^Hl^^^nff. + ^JufpA- (3-7) 

The quality of this approximation is discussed in appendix |^, where the additional corre- 
lations resulting from eq. ( pl6|) are found to be negligible. Note that rule ( p.7|) implies that 
different exciton operators are uncorrelated. Therefore the cross terms appearing in (|3.1|) do 
not contribute to the strength function, which may consequently be obtained by summing 
eq. (|3.3|) over k, a and q. With the approximate correlation rule one easily arrives at: 



'Sl;-^{E', E)=v^ E E (II (A^:)' (Ah_) I 'h?8{E' - El)5{E - E^) . 

{Hc;H+P+;H-P-} {LT;HP} 

(3.8) 

The contribution of the second term on the r.h.s. of eq. ( |3.7| ) is restricted to the diagonal 
elements of the matrix and is therefore neglected relative to the first one. By eq. ( p.8|) 
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the problem factors into two aspects: All information that is specific for the interaction is 
contained in f^. In the present paper this factor is taken for granted. The remaining sum 
is common to all interactions of body rank K. It represents the phase space aspect of the 
problem, the propagation of the i^'-body force into the A-body space. 
The summation over Hq, which appears for all operators with k < K is now trivial and 
yields the factor {^^^ . 

(iii) The sum over {LT; PH} is rewritten in a form which one can call the separation 
of actors and spectators. The configurations {i/+P+} and {iJ_P_} that appear in the 
operators A are called actors. They must appear in {W and | respectively, if the matrix 
element in eq. (|3.8|) is to be different from zero. The remaining configurations that are 
possibly present in (|^| and | ^) are called spectators. As proven in appendix one finds 

E (II(a^:)'(a^:)I^)^ /(AT,i/,p) 

{LT; HP} 

= f{{H~h_,H+),iP-p_,P+),iH-h_,H_),iP-p_,P_)) , (3.9) 

{H-h_.P-p_) 

where / is any function that is completely symmetric in the indices contained in L, T, H and 
P (for each group separately). The sum on the r.h.s. of eq. ( p.9|) runs over the spectators. 
The notation {H — h^, P — p^} means a string oi h — indices of holes and a string of 
p — p- indices of particles that observe the restrictions of eq. ( |2.22 ). In eq. (|3.9| ), there is the 



additional restriction 7^ H_P_H^Pj^ which means that none of the indices in {H — h_,P - 
p_} is allowed to coincide with any one of the indices in H_P_H^P^. By help of eq. ( p.9| 
one can simphfy eq. (|3.8|) as follows 



{H+P+;H.P^}iH-h_^,P-p^} 
jtH^P^H_P_ 



Y: j^y ^(E' - y - e^;)5(E - y - e^-_)5(y - 6^T_ - S), (3.10) 



{H+P+;H_P_} {H-h_,P-p_} 
7^H^P^H_P. 



where we have used E^J'^_ = e^_^_ — S. As already mentioned, the sum over the spectators 
is not independent of the sum over the actors: The spectators are not allowed to occupy the 
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exciton states of the actors. This is only a weak condition if the dimensions of the single 
particle and single hole spaces are much larger than the number of excitons that occur. If 
the energies E and E' are not too high, this will be true and one can treat the sums in 
eq. (|3.1CI| ) as independent. This is called the dilute gas approximation (DGA). It allows us 



to express the strength function by the convolution of partial densities of states 

^l';^iE',E)=v'(J^'^^ J dyp,,,,,{E'-y + S)p,_^^_{E-y + S)p,_,_^^_^_{y) ■ (3.11) 

As discussed in detail in appendix ^ the conditions for the validity of the DGA and for the 
applicability of the approximate correlation rule ( ^.T] ) are essentially the same. 

B. Propagation into the A-body space 



In section [III A| , the matrix element was averaged over configurations with given exciton 
numbers. This result is useful if pre-equuilibrium reactions are studied. In equilibrium CN 
reactions one asks for the average ^ ^ over the full A-body space which is defined as 

^ (E', = < V2 > p {E')p (E) (3. 12) 

= E E m^HmE'-EDSiE-E^) 

t,p=max(A,0) {LT; HP} 



in close analogy with eqs. ( p.lD and ( p.2|) . The summation of the partial strengths up to 
p,t = A is formally correct although in the energy range we are interested in (and committed 
to because of the DGA) by far not all subspaces come into play. The high energy subspaces 
are excluded by the multiplication with (5-functions. As discussed in the last section and 
in appendix different exciton operators are uncorrelated in the framework of the DGA. 
We therefore obtain the strength S{E', E) as a sum of the contributions by the operators 

SiE', ^) = E ''S^''''^ E) , (3. 13) 

kaq 

in obvious notation. These contributions, in turn, are easily expressed by the strength 
functions (|3.3|) : 
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min(y4,A— a) 

E)= Y: KfajE', E) , (3.14) 



p=max(A+fc-42:,2^) 



{p+a,p) \ 

L^ 

2 ' 2 



where the condition t = p + a has been used to evaluate one of the sums over the subspaces. 
The lower limit of the remaining sum guarantees that p > p^, h > (see below eq. (|3.3| )) 
and p> A (see eq. ( p.2C )). The upper limit ensures that p,t < A. We introduce the index 



of summation i = p — and the strength becomes: 



<l+\a\ 



A 2 



fcr,(a,g) 

(i+p+,i+p_) ^ 
i=max(A+fc— (jr,0) 



^S^-^'^\E',E)= Y: ^Zt,+,.)iE',E) . (3.15) 



The second form of eq. ( p.lO| ) then yields 



^^(«'^)(E',£;)=^;2 M Y: E fdyp,,,,,{E'-y + S)p,_^^_{E-y + S) 

{A + k-qfi) ^ ' 

^ 6{y - e\_^_,_^- S) . (3.16) 



Here, the partial level densities of actors are the same as those in eq. (|3.11|) . In eq. (t3.16|) 
the expression 

E E 5(y - e^-A' - 5) . (3.17) 

p=max(A',0) {P-A',P} 

appears with A' = A — ^^^^ and A' = A + k — q = A + z. We want to compare R with the 
total density of states p of the A-body system, the definition ( p.37|) of which is quite similar 



to the expression (|3.171) . The comparison with p{y) is possible if we exploit the invariance 
of the nuclear level density under shifts of the exciton vacuum. We recall eq. ( |2.39| ): 



Piy)= E E S{y ~ e^p_^, - S') . (3.18) 

p=max{A',0) {P-A',P} 



Expressions (|3.18|) and ( p.l7|) differ in two respects 



1. The subspaces with p = A' + 1 . . . A do not appear in eq. ( p.l7|) . 

2. The argument is shifted by S' — S. 
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In the energy range of typical CN reactions the first point is irrelevant so that the only 
remaining difference is the shift of the argument: 

R(?/) = p{y + - + £c'A' - £cA) (3.19) 

In the approximation of equidistant single body levels we find the relation: 

R{y) = p{y--{'-^ + ^S\. (3.20) 



g\ 2 

At first sight one may be surprised to find that this relation depends on A, which a priori we 
may choose arbitrarily. On the other hand, however, the choice of A determines the quality 
of the DGA: it is best if those subspaces that contribute most to the strength are made 
up of as few excitons as possible. The more states are excluded by the Pauli-principle, the 
larger is the error in eq. (|3.16|) . Evaluating the integral, we choose A = 0, which optimizes 
the DGA. Hence, 

^S^'^^^)[E\E) = v'(^^^^ J dyp,_^^^^{E'-y)p,_^^_{E-y)p{y~E,) , (3.21) 

with 

= ^ ' , (3.22) 

which is the energy needed to create z particles (negative z corresponds to the creation of 
holes). 



IV. RESULTS 



A. Transition rates in the exciton model 



Formula (|3.11|) can be applied to the exciton model of preequilibrium nuclear reactions. 



In this model, which was formulated by Griffin |^0[ and has later been refined by several 
authors PT]-|23], an important concept is that of transition rates between subspaces of dif- 
ferent exciton number. These appear in the Master equation for the time dependence of 
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the system. Transitions are assumed to be caused by the residual interaction of a strong 
two-body force V. The rate for going from the subspace Up to Wp+a is given by 

\p^p+a{E) = ^ < > p,+„,p+„(i?) (4.1) 

and is related to the strength calculated in section |111 A| according to 



\^,UE)p,^,{E) = ^-^Y: E) . (4.2) 

Note that the residual interaction consists by definition of the operators with exciton rank 
k = 2 that appear if V is expressed in the exciton picture using the Hartree-Fock single exci- 
ton configurations. Invoking the Ericson densities ( |2.40|) and A = 0, the convolution (|3.11|) 



can easily be evaluated and leads to the transition rate 



Xr^pUE) = fv^ E M:fiE) , (4.3) 



with the density of final states 



-1 



where N = p + h, see above eq. ( ^.1[ ). Because of the last binomial only operators with 
/c > |a| contribute. The final state densities that appear here in the context of propagation 
of the defining GOE matrix elements have been obtained earlier for = 2 by combinatorial 
arguments on the states accessible in two-body collisions [^,^,|l^. Originally, was a fit 



parameter. It is identified here as the average square of the antisymmetric K-hodj matrix 
element. 



B. Spreading width in A-body space 



Formula ( |3.21|) can be evaluated by inserting Ericson's expressions for the partial and 
the CTA for the total density of states, see eqs. ( |2.4U| ) and ( |2.34| ) . The operator A:V(a,g) 
then yields the strength: 
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We have chosen yo = E in eq. ( p.34| ). As mentioned in section [11 D| , below E' 10 MeV 



one should use the temperature T tabulated in ref. [|T6|. For > 10 MeV, the temperature 



should be determined from eq. ( |2.35D , again with i/q = E In experiment, of course, the 



effect of the interaction V as a whole is measured. The spreading width ( ^.1[ ) is obtained by 
summing eq. (|4.5|) over k,a and q after dividing through p{E). This gives: 

Usually a two-body ansatz is made for the interactions between nucleons. We therefore 



explicitly carry out the summation in eq. ( [4.6|) for this case. The factor exp(— 2;(2;+ l)/2gT) 



by which the total spectator density deviates from the nuclear level density — close to unity 
for typical values of T, g and E — is ignored here in order to analyze the general properties 
of the spreading width. One finds: 

T^{E) = 271 v^g^T[2dH + 5{gT)''] , (4.7) 

which is a remarkably simple result. We emphasize two aspects of it: 

(i) According to ref. this is a constant as a function of E below E ^ 10 MeV 



because the nuclear temperature should then be independent of E. At higher energy the 
leading term behaves as E^^'^. Altogether this amounts to quite a weak energy dependence 
of for moderate E — especially if compared to the exponential energy dependence of 
p. This result is in qualitative agreement with the systematics of the spreading widths 
pertaining to isospin violation 0. Eq. ( |4.6| ) shows that the leading energy dependence of 
will be if the body rank of V is i^' instead of two. This demonstrates that the 

reason for the weak energy dependence of the experimental is the two-body character of 
the symmetry breaking interaction: With increasing excitation energy the complexity (in 
terms of excitons) of the states increases. This decreases the fraction of states that can be 
connected by an interaction of low rank, hence the local average square ^ ^ decreases. 
Consequently, the product of ^ ^ and p varies slowly. 
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(ii) The contributions ^T^ of the exciton potential (k=l) and the exciton scattering (k=2) 
to the spreading width are given by the first and second term on the r.h.s. of eq. ( [4.7[ ), i. e. 
by the terms proprotional to T and respectively. Inserting the typical value of T = 0.5 
MeV for E < 10 MeV ^ and g ^ A (13 MeV)"^ as well as dh = A (which optimizes the 



2pi 



DGA), one finds 

1 : ^< 10 MeV 

(4.8) 

lOMeV/E : E>10MeV. 
Thus the rank one and rank two exciton interactions contribute about equally strongly to 
rJ- in the energy range of typical CN reactions. 

C. Discussion 

Within the framework of the statistical model the spreading width due to an arbitrary 
K-hodj force has been calculated for a compound nuclear reaction of an A-hodj nucleus. 
The spreading width is the adequate measure for symmetry nonconservation in complex 
many-body systems 0J2^. It measures the extent to which the states are smeared out due 
to the presence of the symmetry breaking force. It also appears as the damping width of 
a simple configuration — such as an isobaric analog state or a giant resonance — into the 
complex compound nuclear configurations. Definite numbers for a particular interaction 
can be given if eq. ( |4.6|) is complemented by the calculation of the spectral average of the 
interaction in i^-body space. Obvious applications are the electromagnetic and weak forces, 
responsible for the breaking of isospin and parity, respectively. Numerical studies in that 
direction are currently under progress. For the case of isospin breaking a weak dependence 
of the spreading width on energy and mass number is known [Q]. In the present study the 
energy or temperature dependence has been related to the body rank of the underlying 
interaction. Since f ^ depends on A, eqs. ( [4. 61 ) and ( |4.7] ) do not exhibit the complete A- 
dependence of the spreading width, which is expected to be weak. The variation of with 
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temperature can be compared with the results of other authors. Assuming that the spread 
of a one-exciton configuration (a Ip Oh state) is proportional to the square of the excitation 
energy p9| , |30| , Lauritzen et al. have concluded a T^-dependence for the spreading width 



in general. Since the decay of the one-exciton configuration is caused by 2V(i,3), we indeed 
find from eq. ( [4. 41 ) its damping to be proportional to E"^ and eq. ( [4.7] ) indeed indicates 



a T^-dependence of the spreading width in A-body space (for transitions caused by the 
strong interaction only). Studying the coupling of surface modes to single particle motion, 
Esbensen and Bertsch found the "elementary damping" to be proportional to which 
in the framework of Lauritzen et. al corresponds to a T^-dependence of the spreading width. 
There is nothing analogous in the present results since we did not consider collective motion 



in the present paper. De Blasio et. al report the damping of giant resonances to be 
independent of the nuclear temperature. Again the reason is that the statistical damping 
of the present study is different from the damping mechanism considered there. This is also 
true for the study of spreading properties of isobaric analogue and Gamov- Teller resonances 
by Colo et al. |Q. This paper indicates, however, that the statistical damping must be 
taken into account in order to explain the widths of the resonances. 

Frequently the treatment of parity violation in CN reactions is restricted to the k = 1 
part of the weak interaction [p5H 39(|. Eq. (|4.8| ) indicates, however, that the contribution of 



the operators with exciton rank two should be included in the calculation of a local mean 
square matrix element. 

Finally we point out the limitations of our method. The Hartree Fock method and the 
particle hole formalism rely on a basis of product states. Such an independent particle 
model can not describe collective phenomena in nuclei. The present results therefore must 
be modified if applied to reactions that involve collective excitations. A second problem 
are the effects of the symmetries respected by the interaction under study. These effects 
will of course manifest themselves automatically when the spectral average in the 7^— body 
space is calculated. The use of partial level densities for the actors that contain all states 
at a given energy irrespective of further quantum numbers, however, neglects possible local 
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effects of the respected symmetries. We give an example to illustrate this complication. 
Consider the weak potential {k = 1) in a nucleus. This operator connects many-body states 
that differ by only one single-body configuration. Parity violation and conservation of total 
angular momentum demand that the single-body states differ in I but not in j. Since single- 
body states with Al = ±1 and Aj = only exist in different shells, they are separated 
by a relatively large energy interval. Consequently, the present local average suppresses 
these contributions to the spreading width. The two-exciton part of the weak interaction 
on the other hand is not subject to this "local selection rule" because it simultaneously 
changes two single-ecxiton configurations. This fact was first pointed out by Lewenkopf and 
Weidenmiiller [0. They estimated that the k = 2 part of the weak interaction dominates 



the local mean square matrix element. In section [IV B| it was found that potential and 



scattering contribute about equally to the spreading width without taking the local effects 
of the symmetries into account. Sufficiently elaborate expressions for the single-exciton level 
densities of the actors before and after the interaction would only overlap in a small energy 
range, and consequently the convolution of Ph^^p^{E — y)pf^_ ,p_{E — y) with p{y — Ez) would 
result in a smaller contribution of the potential to the strength. 
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APPENDIX A: GOE CORRELATIONS IN THE EXCITON PICTURE 

In this appendix, the differences between the exact ( |3.6| ) and approximate ( |3.7[ ) ensemble 
average are discussed. We find that eq. ( |3.6| ) leads to two types of additional correlation 
that do not appear in its approximation. Considering two examples, it will also be found, 
however, that these additional terms are of the same order of magnitude as those neglected by 
invoking the dilute gas approximation. Since the central formulae of section |III A| and [III B| 
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rely on the DGA, it is not necessary and indeed would be inconsistent to take these quantities 
into account in the evaluation of the strength function. 
The correlation coefficient for two operators is 

Cfcfc' = (II fcV(a,g)| fc'V(a',g')l = (^1) 

{Lc;L_T_;L+T^} 

{HcH^P^\y\HcH^P_){LcL^T^\y\LcL^T_){l\ (a^^)^ (a^:) (aIi) |^) . 

Estimating this expression by its ensemble average, we obtain: 

iLc;L_T_;L+T+} 

Since hole and particle configurations by definition do not intersect, the Kronecker symbol 
implies p+ = t+ and p- = t_. Addition and subtraction of these equations yields q = q' and 
a = a'. In other words: two operators that differ only by their rank are correlated. This 
is the first difference to the ensemble average ( p.7| ) used in section |III A| , were it was found 



that different operators are uncorrelated. The second type of additional correlation we find 
in eq. ( [A2D results from the intricate restriction pattern in the sum over hole configurations: 
The restriction only applies to the indices within the groups Lc and L_. The condition for 
nonzero correlations, however, is that the set LqL^ of lc + I- indices coincide with the set 
HcH_ of he + indices. We adress these two cases one after the other. First, let k' be 
equal to + ^ with > 0. The sum over T+, T_ may now readily be carried out: 

c^u^i = v' E &HfH\{l\ (AS:)' K:) I 'h){1\ (Af:)' (a^:) I S) . (A3) 

{Hc;H-P-;H+P+} 

This expression may be simplified rather easily if we consider the special case of lc = 0. 
Then He = S and we obtain 

Ck,kH = E (II (A^:)' (A^:) I h){1\ (a^^.,)' (a^^_) I 'h) ■ (A4) 

{B;H^P.;H+P+} 

For further simplification of this expression let us introduce the product of hole number 
operators 
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(A5) 



for S = {f3i . . . (3^). This operator has the property 



E(Lh)I^) 

{H} 



(A6) 



where a term containing the factorial of a negative number is understood to be zero. Note 
that 



(A7) 



One therefore obtains 



) E (ii(AS:y(AS-_)ig)^ 

/ {H_P_;//+P+} 



(A8) 



The quantity h — h- is the number of hole spectators, cf. section III A . For the correlation 
coefficient of an operator with itself, we find 



kk = v 



ft) E {Il(AS:)'(AS:)lg)^ 



(A9) 



so that 



\Ck,k+i\ _ 



Ckk {h~h^-OW. 

Under the assumption that K <ti A = dh, which implies ^ <C dh, this yields 

|Cfc,fe+5l ^ / h — h^\^ 
Ckk \ dfi I 



(AlO) 



(All) 



This type of correlation may therefore be neglected if the number of spectators (in hole 
space) is small compared to the dimension dh of the single hole space. This is certainly true 



if /i <^ dh, which is the condition for the DGA introduced in section [III A 



The second type of additional correlation appearing in eq. (|A2|) shall be illustrated using 
the operator iV(o,o) of a two body interaction (cf. Tab. I). We find for its correlation 
coefficient: 
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— 1.2 \^ A^'^ A'^'^ /^l^t^ I P\/T\.i _ I P\ 



MS 



E(Il4a.l I) ((I|aia,| + (I| ^)/. - 2(I|ata,| ^)) . (A12) 
76 



In the brackets, the first term is of the type discussed in sections [III A| and [IIIB| . The second 
term is suppressed relative to the first one by a factor hjdh and by the fact that it contributes 
only to diagonal elements. The third term is suppressed by a factor ^jd^. Obviously, it is 
justified to neglect the second and third term as long as the DGA makes sense. Since out 
of the nine operators of a two-body interaction which actually contribute to the strength 
function iV(o,o) is the only one that shows correlations of this type, eq. (|3.7|) seems to be a 
very good approximation. 



APPENDIX B: PROOF OF EQUATION (Q 

In eq. (|3.8| ), consider the special case of p+ = /i+ = /;,_ = 0: 

-M-E(Sl(A^)l2)V(P). (Bl) 

Let P, P_ be 

P = (n . . . rp) and P_ = (si . . . Sp„) (B2) 

and / a function that is completely symmetric in the arguments t\. . . r^. The restricted 
sum of eq. (|B1| ) can be written as the unrestricted sum 

-M = ^E(S|a..---a.,_at. ...atj$,)2 /(P) . (B3) 
V- p 

The matrix elements vanish unless the si . . . Sp_ all appear in ri . . . Vp. Consider a term that 
satisfies this condition. There are (^^ ^ ways in which the indices r that agree with one of 
Si . . . Sp_ can be distributed over the postions 1 . . .p. Therefore, imposing the requirement 
that the first p_ indices ri . . . rp_ should agree with si . . . Sp_ (up to a permutation) one 
obtains 
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A^ = ^fM E E (S|a.,...a.,_at_...a|,at^_...atj?,)^ /(t,...W,P_) 

P- \P—/ ti...tp_p_ ri...rp_ 

^Si...Sp_ 

p- \P-/ {P-P-) 

Here, the restriction ti . . . tp_p_ ^ si . . . Sp_ means that none of the indices ti ,i = 1 . . -p—p-, 
is allowed to coincide with any one of the indices s^, /c = 1 . . The short hand notation 
{P — p^} ^ means the same. Eq. ( P4|) is obviously the same as 



E illH-'-f fiP-P-,P-)= E ^l-P- fiP-P-,P-): (B5) 

which is easily generalized to eq. 
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FIGURES 

FIG. 1. Two-body interaction in the exciton representation: Diagrams with k = 1,2. 
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TABLES 



TABLE L Two-body interaction in the exciton picture. 



Classification Operator Diagram 
oV(o,o) +E{/3.}(/37|v|/37) Q Q 



iV(-i,i) +E/37«(/37|v|/3s)a^a, 



iV(i,i) +E/3-y.(/3^|v|/37)a!at 



2V(-2,2) + E{/37}{sf}(/37|v|si)a^a/3a,at 



2V(_i,i) + E{/3*}«7(/3^|v|75)4a*a/3a^ 



iV(o,o) -E/3*-y(/3^|v|/37)ata^ Q_ 



iV(o,2) +E/3t.(/3^|v|/3f)atat Q_ 



A" "A 



A 



TABLE I. (Continued) Two-body interaction in the exciton picture. 



Classification Operator Diagram 



2V(_i,i) + Y.i3{su}t{Pt\v\su)alapa, 



2V(o,o) +E{7e}{/3«}(7e|v|/3^)aJaJ- 



2V(o,2) -E7f/3«(75|v|/3f)ataja^at 



2V(o,4) + E{M{««}(5«|v|ii^)aiatata„ 



2V(i,i) + Es{i3-y}s{^s\v\f3j)alala\ag 



2V(i,3) +E„/3{«f}(«*|v|/3w)aJataJ 



a,, 



"A 



V 



V 



2V(2,2) +E{/37}{«f}(^*|v|/37)aIataJat \^ 



